A MECHANICAL MODEL FOR FOURIER'S LAW OF HEAT CONDUCTION. 



by David Ruellef. 



Abstract. Nonequilibrium statistical mechanics close to equilib- 
rium is a physically satisfactory theory centered on the linear 
response formula of Green-Kubo. This formula results from a 
formal first order perturbation calculation without rigorous jus- 
tification. A rigorous derivation of Fourier's law for heat conduc- 
tion from the laws of mechanics remains thus a major unsolved 
problem. In this note we present a deterministic mechanical 
model of a heat-conducting chain with nontrivial interactions, 
where kinetic energy fluctuations at the nodes of the chain are 
removed. In this model the derivation of Fourier's law can pro- 
ceed rigorously. 
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0. Introduction. 

To understand the transport phenomena of nonequilibrium thermodynamics from the 
point of view of microscopic dynamics (say classical mechanics) is a serious challenge. For- 
mally, this is a problem of linear response, solved by the Green-Kubo formula, which is 
basically the result of a first order perturbation calculation. But this perturbation calcu- 
lation is an uncontrolled approximation, as remarked by van Kampen [22] . A fundamental 
derivation of Fourier's law for heat conduction remains thus an open problem, as repeatedly 
pointed out by Lebowitz (see for instance Bonetto et al. [4]). 

Let po be a microcanonical equilibrium state, which is an invariant probability mea- 
sure for the microscopic dynamics of the physical system of interest. The linear response 
problem involves finding the physical state p which replaces po when a small change is made 
to the microscopic dynamics. Problems of this type are mathematically well understood 
when the "microscopic dynamics" corresponds to uniformly hyperbolic smooth dynamics 
on a compact manifold M and p is a so-called SRB state on M{*). Uniform hyperbolicity 
is however too strong a requirement in the physical situation of interest here. In fact, as 
a consequence of the spatial extension of our physical system, there appear a number of 
"central directions" for the dynamics, i.e., we have to deal with partially hyperbolic dy- 
namics. We shall see how this happens in a model discussed below, and how some results 
of Dolgopyat [8] can be applied to a situation where the dynamics is partially hyperbolic, 
and SRB states are replaced by w-Gibbs states. 

To be specific, the purpose of this note is to discuss a deterministic mechanical model 
which exhibits realistic behavior for heat conduction. To obtain our model we start with 
a Hamiltonian chain of -|- 1 nontrivially coupled mechanical systems (nodes) , and we fix 
the temperatures To,T/v of the endpoints of the chain. The kinetic energies of the inter- 
mediate nodes fluctuate. In our model we remove the fluctuations by thermostats that fix 
the intermediate temperatures. We define a stable temperature profile by requiring that 
the intermediate temperatures be such that, for each thermostat, there is no net flux of 
energy in or out of the corresponding node. If this requirement is not satisfied we expect 
the intermediate temperatures to move towards the stable temperature profile when the 
thermostats are removed (see Section 4). Our model allows thus to determine the interme- 
diate temperatures for a heat-conducting chain. Removing the kinetic energy fiuctuations 
of the nodes of the original Hamiltonian chain is an uncontrolled approximation, but once 
it is accepted one can proceed rigorously. One of our results will be Fourier's Law: the 
amount of energy transported by the chain is asymptotically ~ N~^{Tn —Tq) for large N 
and small Tjv — Tq (see Section 5 for a more prudent statement). 

(*) There is a large literature on this subject. See for instance Young [23] for a general 
discussion of SRB (Sinai, Ruelle, Bowen) states. Linear response for uniformly hyperbolic 
diffeomorphisms was established by Katok ct al. [16], sec also Ruelle [18]. The flow case, 
which is most relevant for physics, is discussed in Ruelle [22], Butterley and Liverani [5] (For 
applications to physics, see Gallavotti and Cohen [13], Ruelle [19]). Further discussion and 
references can be found in Dolgopyat [8] who discusses extensions to partially hyperbolic 
dynamics. An idea of the variety of nonhyperbolic dynamics is given by Bonatti et al. [3]. 
See also in this respect Baladi and Smania [2], and Ruelle [21]. 
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The treatment given here strives at conceptual clarity rather than generality. At the 
cost of increased mathematical complexity (or new ideas) one could probably deal with 
much more general situations than the simple model discussed below. Note that there 
are a number of rigorous papers related to the statistical mechanics of heat conduction, 
but using approaches different from that presented here. This includes work by Eckmann, 
Gallavotti, Hairer, Jaksic, Liverani, Pillet, Rey-Bellet, Young, etc. (see in particular [9], 
[1], [10], [11]). See also a study with stochastic thermostats [4a], [4b], and a promising 
investigation by Dolgopyat and Liverani [8a] of the macroscopic behavior of a coupled 
lattice of strongly chaotic microscopic subsystems. 

Acknowledgments. 
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paper [8]. I am also indebted for a number of discussions on linear response in smooth dy- 
namics and nonequilibrium statistical mechanics to Viviane Baladi, Jean-Pierre Eckmann, 
Giovanni Gallavotti, Vojkan Jaksic, Joel Lebowitz, and Lai-Sang Young. 

1. Our model: the time evolution (/*). 

Our model is a chain of N + 1 nodes. The nodes, before coupling, are assumed to 
be Hamiltonian systems described by geodesic flows on a compact n-dimensional Riemann 
manifold M, with n > 2. The Hamiltonian Hj of the j-th node is thus the kinetic energy. 
In local coordinates 

where we have written pj = (pju), Qj = {qj), = Yliuv d"""" {.^j)<^ubv, and (fif'™) is the 

inverse of the matrix {quv) defining the metric (the nodes have mass m). This choice of 
Hamiltonian system gives examples where the time evolution is an exponentially mixing 
Anosov flow. Furthermore there will be a simple relation between the kinetic energy and 
the temperature (see Section 4). 

We write Xj = (pj, q^) and let (/j) be the geodesic flow (xj, t) i— )■ /jxj restricted to the 
energy shell Sj = {xj : ii"j(xj) = Kj} with Kj > 0. Deflning S = xf^^Sj C T*{M^+^), 

we also let /* = x^q/* on S. The time evolution defined by the Hamiltonian Ylo ^j, 
when restricted to S", is thus (/x)- 

We introduce now a coupling between the node j and its neighbors, given by a force 
XX j G T*, M where A G R. We require that there is a smooth function W : M x M ^ Ti 
[satisfying condition (1.3) below] such that 

x,=x-+x+ 

where Xq = X^ = 0, and the other X^ are given by 

X- = -d^Wici^_uCLj) , X+ = -aq^.W^(q,-,q,+i) (1.1) 
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[For simplicity we do not introduce a self- force X° depending only on qj.] The time 
evolution corresponding to the coupled Hamiltonian 



N , N 



2 m 

3=0 j=l 



is given by 



Write now 



± fPju \ ^ f -Zvwi99'''"/dqJ)p,.p,^/2m + XXj^ 



OLj = {Xj,pj) = CK- +(5;+ where = {X^,pj) 



aj = aj/ (pj, pj) = a- + a+ where a- = 6l- /(pj, Pj) 

We specify our model to correspond to the following coupled time evolution (/*) on S for 
the A/" + 1 nodes: 



(1.2) 



The choice of the aj is such that for the coupled time evolution we have 

uv uv 

for j' = 0, . . . , A^. The term —Xajpju in (1.2) is called an isokinetic thermostat (introduced 
by Evans and Hoover, see [12], [15]): it keeps the kinetic energy of node j fixed to a value 
Kj for each j. [This is a physically reasonable thermostat, especially when n is large.] 
Therefore (/*) is a time evolution on S as announced. (We shall complete the description 
of our model by making a specific choice of the kinetic energies Kj, see the definition of a 
stable temperature profile in Section 5.) 

Let us define 

Pjidy^j) = dpjdcij 

where dpj is the normalized volume on the sphere {pj : Yluv 9^^i^3)PjuPjv = 'imKj}, and 
dqj is the normalized Riemann volume on M; pj is thus an an ergodic measure for (/j) on 

Sj. We also write x = (xq, . . . , xjv), and define px (dx) = YljLo Pj{d:x.j). We assume that 

J dcij-iW{cij-i,qj) = j dc[jW{cij-i,qj) = (1.3) 

Note that, since dpj is invariant under pj — Pj, we have 

Px(«,^) = (1.4) 
4 



Furthermore, 

Px o /x)) = unless al is a+ or aj_^^ (1.5) 

Px{{c(J){c(k° fx)) = ^ unless q;| is a" or Q!t_^ (1.6) 

To see this, note that /* does not mix different x^, that px is (/* )*-invariant, that 
J dpi{X^,pi) — (invariance of dpg under pi i— )■ — p^), and that J dqi±i{X^,pi) — 
(because of (1.1) and (1.3)). 

2. A perturbation argument. 

Since pj has a smooth density on Sj, we see that it is an SRB state for (/|), and also 
that Px is an SRB state for (/* ) on S. We shall now assume that the (/j) are exponentially 
mixing Anosov Rows. (Since (/j) is the geodesic flow on a compact Riemann manifold M, 
this is the case if M is a surface of negative curvature, see [6], [7].) We refer the reader 
to Appendix A for a precise definition and a proof that (px)(/x)) is also exponentially 
mixing. In particular {px, (fx)) is ergodic. 

To study the physics corresponding to the perturbed time evolution (/*) defined by 
(1.2), the existence of an SRB state p would be desirable. But, since (/*) is not uniformly 
hyperbolic, such a state need not exist when X ^ 0. It is however possible, following an 
idea of Dolgopyat [8], to study a perturbation of px to w-Gibbs states, as we now explain. 

A dynamical system (/*) on the compact Riemann manifold S is said to be partially 
hyperbolic if there is a continuous invariant splitting 

such that, for suitable constants C > 0, 6*1, ^2 < 1, and all t > 0, we have 

(VveE'^) \\{Tf-*)v\\<C9{\\v\\ 

(Vx e S) \\{Tf)\E^M{Tr')\E^\\-' < Ce\ 

One can then define local unstable manifolds W^; the corresponding global manifolds form 
a continuous foliation W" of S with smooth leaves, which is tangent to E'^ (see [14]). 

An (/*)-invariant probability measure p on 5" is called (by Pesin and Sinai [17]) a 
u-Gibbs state if the conditional measures on the local unstable manifolds have a density of 
a certain canonical form. The tt-Gibbs states are precisely the (/*)-invariant probability 
measures p on 5 which are absolutely continuous with respect to the foliation (i.e., if 
XnWx has leaf Lebesgue measure for each local unstable manifold Wx, then p{X) = 0). 
See Dolgopyat [8] for a discussion of w-Gibbs states from this point of view. If p is SRB, 
then p is also w-Gibbs. If £ is a probability measure absolutely continuous with respect to 
the Riemann volume on S, and 
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then any weak limit of It when T — > oo is a tt-Gibbs state. In particular, the set of tt-Gibbs 
states is nonempty. 

A tt-Gibbs state, while (/*)-invariant, may have a natural decomposition into states 
corresponding to slow time-oscillations of our physical system. This possibility makes u- 
Gibbs states more flexible objects than SRB states for the description of nonequilibrium 
steady states. 

In the physical situation that we want to discuss, (/* ) is partially hyperbolic because 
the (/j) are Anosov flows. The perturbed time evolution (/*) is thus also partially hyper- 
bolic, provided A is small enough (see [17]). For A 7^ there need not be an SRB state, 
but limits of It are tt-Gibbs states, and the following applies: 

3. Proposition. 

For A sufRciently small, if p is a u-Gibbs state for (/*), and A a smooth function on 
S, we have 

p{A)-p^{A)^X{n-l) dTpx(E«,)(Ao/J)) + o(A) (3.1) 



This result is a corollary of a theorem of Dolgopyat [8] , which applies to the time 1 map 
fl- because of our assumptions. Specifically, is a rapidly mixing (in fact exponentially 
mixing) Anosov clement in an abelian Anosov action on S. Therefore, for sufficiently small 
A, if p is a u-Gibbs state for (/*) (hence /^), Theorem 1 of [8] yields 

piA)-p^iA) = XuiA) + oiX) (3.2) 

where A <^{A) is a linear functional on smooth functions. We refer to [8] for the 
definition of abelian Anosov action and other details including weaker conditions on A. A 
proof that is rapidly mixing in the sense of [8] is given in Appendix A. 

Dolgopyat in [8] gives an explicit expression for the linear functional uj. Since 
preserves the volume element pxidx) on S, Dolgopyat's expression can be simplified to 

00 

uj{A) = -Y,Pxi[d^yY][Aof^]) (3.3) 

fc=0 

where the divergence is taken with respect to px, and 

We refer to Appendix B for the proof of (3.3). Finally, we can compute divY in our case; 
this is done in Appendix C and yields (3.1). [] 

4. Energy transfers. 
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Taking A = in (3.1), and using (1.4 — 1.6), we have 

nOO 

p{af) = X{n - 1) / dT{pjPj±i)[{af + af^,){af o /J)] + o(A) 

^0 



Since on Sj we have = 2mKja^, we obtain 



piaf) = - / dripjpji.maf + pj±,af^,)iaf o /J)] + o(A) (4.1) 

'IT' Jo 

where j3~^ = 2Kj/{n — 1) is defined to be the temperature associated with the kinetic 
energy Kj. 

Using the invariance of pj under (p, q) i-^ (— p, q) we have 

/ dT{pjPj_i)[{a+_-^){aJ o /J)] = / dT{pjPj^i)[{a+_-^){aJ o f'^)] 
Jo Jo 



f 

J —c< 



dT(pjPj-i)[(a+)(a. ofl] 



and similarly with interchange of + and — , so that (4.1) becomes 

A f°° 

P^^^^ = ^ cZT(p,p,±i)[(/3,a± + /3,±iaJ±i)(a± o /J )] + o(A) (4.2) 

Let ((/)*) be the geodesic flow with unit velocity on M, and the corresponding 
absolutely continuous invariant measure. Writing VF(xi , X2) instead of VF(qi , q2) we define 

*(Pl5P2,Ti - (71, T2 - 0-2) 

= j p^{d^i)p4d^2)w{r''''^i,r''''^2)w{r'^'xi,r'^'^2) 

^ijipi,P2,Tl - (71, T2 - (72) = daidrj'^(pi,P2,Tl - (7i,T2 - (72) 

and 

/+00 
dT^ij{pi,P2,T, t) 
-oo 



Then wc may write $ = $11 = —$12 = —$21 = ^22 ^ 0. Because of the exponential 
decay of correlations, $(^1,^2) depends smoothly on pi,p2 > 0. Note also that 

Hpi,P2) = Hp2,Pi) , $(p,p) =p$(l,l) 
It is easy to check that 

/+00 
dr {pj-iPj)[{a+_^){a+_^ o /J )] = $(|Pj-i|, \Pj\) 
-00 



/. 
/_ 



+ 00 

+ 



dr{Pj-iPj)[m_i){a. o/J)] = -$(|p,_i|,|p,|) 



dr (Pi-ipj)[(a,- o /J)] = -$(|p,-i|, |pj|) 

c/t(pj_iPj-)[(«7)((«7 °/x)] = *(|Pi-iUp7) 

-oo 

In particular, (4.2) gives 

By definition of o:^, the average total transfer of energy per unit time from the node 
J — 1 to the node j is thus 



72 



-[p{a-) - p{a\_,)] = -oiPj - ^,-i)$(|p,-i|, |p7) + o{X') (4.3) 



m ■' ■' m 

ior j = 1, . . . , N. Note in this respect that 

P(Q!/_1 + Oij) =0 

This is because 

+ aj = (X+_i,pj_i) + (Xr,pj.) 

= (a,,+a,JW^(q,_i(ti),q,(t2))|,^=,^^, = Wcb-i(t),q,(t)) 

and the average of the right-hand side over the (/*)-invariant measure p vanishes. 

The average transfer of energy per unit time from the thermostatting force —X{aJ + 
Q!^)pj/m to the node j can be obtained by a similar calculation. It is 

——p(aj + af) 
m ■' ■' 

= |Pi-il) + (^,- -^,-+i)*(|p7, |p,-+i|)] +c,(A^) (4.4) 

This is also minus the initial rate of heating of the node j in the absence of a thermostat. 
5. Stable and approximate stable temperature profiles. 

It is natural to define a stable temperature profile (STP) by fixing f3o and /^at, and 
requiring . . . , I3n-i to be such that pi&J + = for j = 1, . . . , — 1. [Equivalently, 
Kq, Kn are fixed, and Ki, . . . , -ftTyv-i are such that p{a~ + a^) = 0.] In an STP, there is 
thus no net energy contribution from the thermostats to the nodes 1, . . . , AT — 1, but the 
thermostats remove the kinetic energy fluctuations at these nodes. 
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We expect that, for given /3o,/3jv, some STP exists, at least when A is sufficiently- 
small, but we do not have a proof of that fact. Technically, what is lacking for a rigorous 
discussion of STP's is a proof of uniformity of o(A), in (3.1) or (4.1), with respect to 
/5o, . . . , (3n in a compact interval* 

This being the case, we shall content ourselves with a discussion of approximate STP's 
defined by 

(/3, - /3,_i)$(|p,_i|, |p,|) = (/3,+i - PjM\Pj\, |p,+i|) (5.1) 

for J = 1, . . . , N — 1. For an approximate STP, (4.4) says that the average energy transfer 
per unit time from the thermostatting forces is o(A^), which is very small for small A. The 
total energy transfer per unit time through our chain, i.e., p{a'^) — p{aQ) is thus, according 
to (4.3) 

= S(/3.-/3.-i)$(|p,--i|,|p,-|) + o(A^) 

I ft 

fOTj = l,...,N. 

We take now for definiteness /3o < /3n- We know that $(pi,p2) > 0, and we shall 
assume that 

max 

when pi,P2 belong to some compact set where the |pj| are allowed to vary: the condition 
^min > expresses that the nodes of our chain are actually interacting. We may then 
rewrite (5.1) as 

where $~ = $(|pj_i|, \pj\), = $(|pj|, |Pj+i|)- We have thus Po < Pi < . . . < Pn, and 
also 

{pj - ^,_i)$„.i„ < ipj - ^,_i)$(|p,_i|, |p,-|) 
1 ^ 1 

= ^T.(^J-^j-^)^i\Pj-ll\Pj\) ^ ^(/5iV -^0)^max 



Therefore 



0<(3j- <^{Pn- /3o)$max/$r 



tends to zero when N ^ oo. 

If K,{N) — /w?){l3j — /3j_i)$(|pj_i|, |pj|) is the total energy transfer per unit time 
through the chain (up to o(A^)) we have thus for large N 



A2 _ A 



2 



n{N) ^ — (/3, - /3,-i)$(|p,|, |p,|) = -^(/3, - /3,-i)|p,|$(l, 1) 

lib fib 



* This uniformity can probably be proved, according to Dmitry Dolgopyat (private 
communication) . 
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hence 

A2$(l,l)|p,| A2$(l,l)^/;^^(^ 
2^1/2 ~ 2A2$(l,l)Vm(n-l) 



so that 



FinaUy 



AT m 

conforms to Fourier's law: 



1/2 _ 1/2 ^ 



iV 

when is large and /3n — f^o small. 



Appendix A: exponential mixing and rapid mixing. 

Let (/j) be a smooth flow on the compact manifold Sj. We say that (/j) is exponen- 
tially mixing with respect to the invariant state pj if for some p > there are 7, C > 0, 
such that 

\p,{{A,ofj)B,)-pM3)Pj{B)j\<C\\A,\\,\\B,\\,^^^^ 

when Aj,Bj e CP{Sj). 

We shall now see that if this holds for j = 0, . . . , A?", then the flow (/* ) on S = Xq Sj 
defined by (xq, . . . , x^r) = (/qXo, . . . , /^xjv) is also exponentially mixing with respect 
to px — Pj- Indeed we shall prove that 

\px{{Ao fi)B)- px{A)px{B)\<{N +l)C\\AU\B\\,e-^' (Al) 

provided A, Be Cp{S). 

Define, for j = 0, . . . , TV + 1, 

Aj{xj,Xj+i, xjv) = J po{dxo) ■ ■ ■ pj-i{dxj-i)A{xo, . ,xn) 

and similarly for Bj. By assumption, for all . . . , x'^y, x^'^.^^, . . . , x^, we have 

I J pj{dx.j)Aj{fjXj,x.'j_^-^^, . . . , x^)Sj(xj, x^'^.j^, . . . , x^) 

-(i,-+i(x;.+i, . . . ,x'^)S,-+i(x;Vi, • • • , x^)| < CWAWpWBWpe--^' 
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hence 



I J Pjid^j)--- J Piv(c^xjv)ij(/]xj,...,/^xjv)-Bj(xj,...,Xiv) 
J Pj{d^j+i) ■J Piv(c^Xiv)ij+i(/|+iXj+i, . . . , f^XN)Bj+i{xj+i, ...,xn)\ 



< C\\A\\p\\B\\pe-^* 



Since 



J po{dxo)---J piv(<ixiv)Ao(/oXo,...,/^Xiv)So(xo, ...,Xiv) = Px((^o/*)S) 
An+i=Px{A) , Bn+i=Px{B) 

we obtain (^.1). 

The case of interest to us is when (/j) is an exponentiaUy mixing Anosov flow with 
respect to the volume pj on 5'^. Then (/* ) is exponentiaUy mixing with respect to the 
volume px on S, and therefore is rapidly mixing in the sense of [8]. This means 
(roughly) that if £u,<t is a probability measure with ry-Holder density cr on a local unstable 
manifold >V", and if A has derivatives in the center direction which are r^-Holder in S, 
then iu^cr{Ao f^) tends to px (A) faster than any n~'^ (with k > 0) when n — >■ oo. To check 
that exponential mixing implies rapid mixing one can approximate 5-measures on unstable 
disks (used to define iu,^) by smooth functions. [I am indebted to Dmitry Dolgopyat for 
explaining this to me.] In fact, spreading the mass of £u,<t by a small distance ~ r along 
stable manifolds, then along center manifolds, we get a probability measure with e- 
Holder density on M such that |4,a(-4 o f^) - £^{A o /^)| < dr"^. The exponent e is 
determined from the Holder exponent of a and of the stable foliations, and we choose e < ry. 
Smoothing (/) G C^to^e we have 1 | < C2r% so that |£^(Ao/J)-£^(Ao/^)| < C^r^ 

Smoothing A to i G we have ||A-i||o < Cgr^, so that |£^(Ao/^) -£^(io/^)| < C^r^. 

Assuming px{,A) = we may also assume Px{,A) — 0, and exponential mixing gives 

K^(^°/x)l = |Px(io/-)0)| < C4r-^e-^- 

and thus 

° /x)l < {C^ + + C'^)r^ + C^r-^e--^^ 

Taking r such that 



we have 

K„,,(^o/-)|<C5exp(-7-^n) 

for large n. In particular, is rapidly mixing. 
Appendix B: proof of (3.3). 
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The proof of theorem 1 in [8] involves the invariant sphtting 



TS = E'^^E^^® E'^ 

associated with (/* ). Here E"^ = ®oE^, E% = ®^E^, Et, = ®oE^-, where E]^ is the 
unstable vector bundle for (/*), the stable vector bundle, and the one-dimensional 
bundle in the direction of the flow, so that [E'^)-^. is spanned by 

e.(x,) = |/;x,/|||/jx,|| 

The bundle E^. is Holder continuous, but in general not smooth. A smooth bundle E^^*, 
C°-close to E?^ , is introduced. The components Z", Z^, Z°-^ of a vector Z will be taken with 
respect to the splitting TS = E\®E%^ ®E'^'. We shall use bundle maps T", T"^ : TS ^ 
TS such that (T"x))z = {{T^p)zr, {T%^))Z = {{T^p)ZY, (r«^(x))Z = ((^x/l)Z)«^ 

Let now 

and define the vector field V and the functions aj on S by 



F(x) = 5](T-rr-(/--x) 



n=0 

j=o 

Then the functional u such that (3.2) holds is given by Proposition 2.6 of [8]: 

N oo oo 

oj{A) = p^idvA) + ((«^- ° /x ")^e,^) - Px (([div^l^^ + T'^F)] o f--)A) 

j=0 n=0 n=0 

where div" is the divergence with respect to the canonical density on W^. 
We claim that we can transform the above formula to 

oo 

oj{A) = -Y,P><{[diyY][Aof^]) (3.3) 

fc=0 

To show that the right-hand sides are equal, we may replace px by its conditional measure 
Pu on a local unstable manifold W", i.e., integrate on W" with respect to the canonical 
volume element. The divergence div is with respect to the volume element px{d:x.), and 
along an unstable manifold it can be naturally factorized in volume elements along the 
as, c, and u directions, so that we have 

-Px ([divy][^ o /^]) = (y- . diA o /^)) - ([div'^-y-if^ o /^]) 
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where 



= ((T-y-) • o ')) - Px ([div'='^(T-y-)] o ']) = • • • 

= Px(((T"^)''y"*) • ^yl) - ^Px([div'="(T'^"(T'**)^-^y''"][ylo /^-i]) 

Using (6) and Lemma B.l of [8] we see that the sum 

|^Px([div^"(T-(T-)^-iy-][Ao/^"i]) 

converges absolutely, so that 

oo oo 

-J2px ([divy] [Aof^])^p^{V-dA)-J2 Px ([div'^"(l^'^" + T^F] [Aof^]) 

k=0 k=0 

and the right-hand side is equal to the expression for oj{A) of [8] reproduced above, 
have thus proved (3.3). 

Appendix C: proof of (3.1). 

In our case 

divy(x) = Adiv /" dt{Tf-t^f)X ^ X [ dt{divX){f-^x) 
Jo Jo 

where X e TS has components 

V 

Therefore 

oo oo „x 

a;(yl) = -^Px([divy][^o/^]) = -A^ / dt p^{[dwX][Ao f'^+']) 

fe=0 fc=0 "^0 

/•oo 

= -A / dtpx([divX][Ao/^]) 
One proves readily the formula dp^ ■ {p.jPj) = {n — l)aj, from which we obtain 



-divX = (n- l)^a 

j^O 



N 

3 



hence 

AT 



/"OO J ' 

{A) = X{n-l) dtp^{C£^,){Aofl^)) 



j=0 
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i.e., we have proved (3.1). [] 
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